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1. Introduction 

It has been known that quasilinear parabolic equations or non-linear reaction-diffusion 
systems arise in physics, chemistry, biology, and other applied sciences. The following three 
equations are the examples of this type of partial differential equations (PDEs). First one is 
called the Korteweg-de Vries equation 


( 1 . 1 ) 


Ut + cUU x + Uxxx - 0 , 


and first encountered in the study of waters, Korteweg [1], denoted by KdV. The other is 
called the Complex Modified Korteweg-de Vries-I equation (CMKdV-I) 



( 1 . 2 ) 
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which arises both in the asymptotic investigation of electrostatic waves in a magnetized 
plasma and in the asymptotic investigation of one-dimensional plane-wave propagation in a 
micropolar medium, Erbay [2]. The last one is the Complex Modified Korteweg-de Vries-II 
equation (CMKdV-II) 


U t -6\U\ 2 U x + U xxx =0, 


(1.3) 


which is another example for quasilinear parabolic equations or non-linear reaction-diffusion 
systems, Ablowitz [3]. Equation (1.2) does not hold the Painleve property but the (1.1) and 
(1.3) do, Mohammad [4]. The equations which have Painleve property may be solved by the 
method of Inverse Scattering Transformations (1ST) and hence they are completely integrable 
[5, 6]. Sometimes it is not easy to solve 1ST problems [3], such as for CMKdV-II equation. 
Therefore the need for an easy and useful method which has to give soliton solutions for a 
given PDE is emerged. An important method is developed by Hirota for finding N-soliton 
solutions of non-linear PDE [5, 6]. 

In this paper, the Hirota's method is applied to the CMKdV-II equation. The Hirota's 
method generally requires the transformation of PDE into homogeneous bilinear forms 
of degree two. Only specific PDEs can be transformed in this way. This means, when a 
bilinear (form) equation can be solved, then N-para meter solution can be obtained as a 
series which self-truncates at finite length. These expansions that self-truncate in this way 
give automatically exact solutions. Self-truncation, however, does not occur for all bilinear 
equations; if it does, then the equation in question possesses multiple soliton solutions. The 
reason for this situation has never been adequately explained. In other words, self-truncation 
which is equivalent to complete integrability would require a connection with the conserved 
quantities of the original equation. 

In this study, it is proven that the CMKdV-II equation has self-truncated Hirota 
expansions. It is shown that there is a direct equivalence between the Af-soliton solutions 
of Hirota's bilinear form of CMKdV-II and the Backlund transformations proposed by Weiss, 
Tabor, and Carnevale [7,8]. 

Now, the question here is where the soliton comes from. Firstly, J. S. Russel in 1834 
recorded his observations of great solitary wave as a mean of developing the mathematical 
properties of a large class of solvable non-linear evolution equations. Solitary waves, solitons, 
Backlund transformations, conserved quantities and integrable evolutions which can be also 
named as completely integrable Hamiltonian systems are in the class of solvable non-linear 
evolution equations. The description of John Scott Russel has aroused among mathematicians 
and physicists one hundred and forty years later, Zabusky [9]. In their paper, they were the 
first ones who defined the solution for the following KdV equation: 


lit + 6UU X + il xxx - 0. 


(1.4) 


For the CMKdV-II Equation, the Hirota's bilinear form is given in Section 2 and the analytical 
one- and two-soliton solutions are presented in Section 3. The numerical procedure and 
results for one-soliton solution are outlined in Section 4. 
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2. Hirota's Bilinear Form of the CMKdV-II Equation 

It is known that the equation 


3 


Li f - 6\U\ 2 U x + U xxx = 0 


( 2 . 1 ) 


is the complex modified Korteweg de Vries II equation (CMKdV-II). Let g and f be the 
complex and real valued functions, respectively, satisfying 


^=- r lsl z = -//»+fl 


( 2 . 2 ) 


By using the transformation above, CMKdV-II becomes 


f gt ftg + 3(fxxgx fxgxx) + fgxxx gfxxx — 0 - 


(2.3) 


Let 


D = Dt = — 
x dx dx' 1 dt dt' 


(2.4) 


then (2.3) becomes 


\d 3 x + D t ]g(x,t)f(x',t‘) \ x , = xt , = t =0 ^(d 3 x + D t )gf = 0, (2.5) 


which is a homogeneous bilinear form and called Hirota's form of CMKdV-II equation. Let 


u= L ^- + u 1 


( 2 . 6 ) 


be the truncated solution of CMKdV-II. Then 


\U 0 \ 2 = ® 2 , 


UqUi + UiUq = -<£>xx, 

6U() X \U 11 - - 6<£>x X U\x= H oxxx + Hot, 

6<Ibiro|iri| 2 - 6<Z> 2 x U 1x = 3(® x U 0xx - G>xxlIox) + Uo(® X xx + Ot), 


U lt -6\U 1 \ 2 U lx + U 1X xx = 0. 


(2.7) 
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Here 11 and U\ are both separated solutions of the CMKdV-II equation. Hence it is an onto 
Backlund transformation of CMKdV-II equation. The system of these equations are called the 
Painlevd relations. There is a relation between the soliton of the bilinear (2.5) and the function 
® of the Painlev6 relations (2.7). Consider the solitons 


having the properties 



U i = 


8 


(i) 


/ (1) 


l^l 2 = (log /)„, 

|LJ|2 = 'S l0§0+l ' Jl12 - 


It can be shown that 


( 2 . 8 ) 


(2.9) 


/ = ®/ (1) , g = H 0 / (1) +®g (1) . (2.10) 

Theorem 2.1. If and g (n> satisfy (2.5) for all n, with 

I g {n) \=f W f^-(fx n) ) 2 , (2.11) 


and if 


/ ( ” ) =0(„_i)/ ( ”- 1) . 


( 2 . 12 ) 


then the resulting equations in ® n _i,LZo and U( n are satisfied by the PainlevS relations (2.7). 
Furthermore 

f (n) = ft®,, (2.i3) 

i=0 

with 

/ (0) = 1. (2.14) 


Proof. When substituting (2.12) into (2.5) and using the Painleve relations when necessary 
yields the claim of the theorem. Using (2.12) succesively we obtain the relation (2.13) that 
completes the proof of the theorem [4]. O 
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By using the usual perturbation method, /V-para meter exact solitary wave solutions of (7) 
can be obtained, Nayfeh [10]. The power series of g and f which are given in Hirota [11] in a 
small parameter e are: 


g = eg + e 3 g^ + e 5 g (5> + ■ ■ 
f = l + £ 2 /( 2 ) +£ 4/(4) + ... / 


(3-1) 


where g and f are the solutions of the (2.5). Then considering the increasing powers of e from 

(2.5) it is clear that 


gxxx + ^f UI = 0. 

(3.2) 

gS,+g, (3| = -(Dt + D i)f m g m - 

(3.3) 

gSx + g! s> = -(d 3 , + d) [/<V> + /‘V 3 ’]. 

(3.4) 

& * = -(Dt + D.) [/<V> + f'V' + /‘V 5 >] 

(3.5) 


and (2.2) yields 

(3.6) 

(3.7) 

f m fS-f m f£. (3.8) 


fS = -*< V, 

fS = -g ( V - g m ’g m - f?f? - f m fS, 
fS - -sr< V - s '"’?® 1 - s“’s ,5> * - 2/P/b - 


where * stands for the complex conjugate. 


3.1. One-Soliton Solution for the CMKdV-II Equation 

To obtain one-soliton solution of the CMKdV-II equation, let's take 


gd) = e® +rF ; 0 + r¥=a l x + b 1 t + c 1 . 


(3.9) 
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Hence 


= re(ai) + I ■ im(ai), 
b\ = re(bi) + I ■ im(i>i), 

Ci = re(ci) + I • im(ci), 

re(t>i) = -(re(ai)) 3 + 3 * re(ai) - (im(fli)) 2 , (3.10) 

im(bi) = (re(fli)) 3 - 3 * im(a!) - (re^)) 2 , 
p 2® 

2 ) = e 

4(re(a 1 )) 2 ' 

and = 0, g( 2n_1 l = 0 for all n > 2 where I stands for the complex number i. Therefore 

f(t,x) = 1 + f 2 e 2<P / g(x,t) = £ c 0+i¥ (3.11) 

4(re(fli)) 


is a solution of Hirota bilinear form (2.5) and 


g(x, t) 


f(x,t) 


(3.12) 


is the corresponding one-soliton solution of the CMKdV-II equation. 

3 . 2 . Two-Soliton Solution for the CMKdV-II Equation 

To find two-soliton solution, let's take 

g(i) = e ®i+m + e 


where 


®i[f,x]=re(fli) * x + re(bi) * f + re(ci) 

¥1 [t, x] = im(ai) * x + im(bi) * t + im(ci) 
<I >2 [t,x] = re(« 2 ) * x + re(b 2 ) * f + re(c 2 ) 
V 2 [t,x] = im(fl2) * x + im(i>2) *t + im(c2) 


(3.14) 
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re(bi) = -(re(fli)) 3 + 3 * re(«i) - (im(ai)) 2 
im(h) = (re(ai)) 3 - 3 * im(ai) - (re^)) 2 
re(i> 2 ) = -(re(fl 2 )) 3 + 3 * re(fl 2 ) - (im(a 2 )) 2 
im(b 2 ) = (re(a 2 )) 3 - 3 * im(a 2 ) - (re(a 2 )) 2 - 

Hence from (3.8), it can immediately be found that 

2 ®! 2® 2 
f< 2 > = —_+ —_. 

4(re(ai)) 2 4(re(a 2 )) 2 


(3.15) 


(3.16) 


Since the right hand side of (3.3) is not zero, a special solution for g (0,> can be found by the 
method of undetermined coefficients. Considering (3.7), is a nonzero real function. The 
right hand side of (3.4) is not zero, and we can take g ^ = 0. From (3.8), it is seen that is 
zero and that = 0 is taken for convenience. Further computations have shown that 

f(2n) = Q' g ( 2 n-l) = Q, „ > 3 . (3.17) 


Hence 


f[t,x] = 1 - 


,2®! [t,x] _ 


„2<b 2 [t,x] 


4(re(ai)) 4(re(a 2 )) 

^(im(fli) - im(a 2 )) 2 + (re(ai) - re(a 2 )) 2 ) 
16(re(a 1 )) 2 (re(a 2 )) 2 [(im(ai) - im(a 2 )) 2 + (re(«i) - re(a 2 )) 2 

_^_ _®1 [t,x] +®2 f t,x] + [t,x] -/■ 1 l f 2 [t,x] 


*2(01 [t,*] +®2[U]) 


[re(aj) + I ■ (im(aj) - im(a 2 )) + re(a 2 )]‘ 
1 

[re(fli) + I ■ (im(ai) - im(a 2 )) + re(fl 2)] 2 

g[t,x] = e ^i[t,x]+^2lt,x]+I^i[i,x]+I-V 2 [t,x] 


p&l [t,x] +O 2 [t,x] -I ^l [t,x] +1 -^2 [t/X] 


{-I ■ im(ai) - re(fli) + I ■ im(a 2 ) + re(a 2 )) 2 0 l [ t , x ] + 2® 2 [t,x] + J’F 1 [t,x] 

ry ^ 

[2re(fl2)(im(fli) - I ■ re(ai) - im(a 2 ) + I • re(«2))] 

{-I ■ im^) - re(ai) + I ■ im(a 2 ) + re(a 2 )) 2 2 ®![f,*] + ®2[u]+PF 2 [f,*] 

[2re(fl2)(im(fli) - I ■ re(ai) - im(a 2 ) + I • re(fl 2 ))] 


(3.18) 
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Table 1: The difference of two consecutive solution values (N = 256, -128 < x < 128, M = 500, 0 < t < 5). 


Ax 

At 

Li 

Lee 

1.0 

0.01 

0.1410966549 

0.6322912259 

1.0 

0.01 

0.0033087030 

0.0105278461 

1.0 

0.01 

0.0000823911 

0.0002140195 

1.0 

0.01 

0.0000012444 

0.0000028444 


is a solution of the Hirota's bilinear equation (2.5) and 


g(x, t) 


f(x, t) 


(3.19) 


is the corresponding two-soliton solution of the CMKdV-II (1.3). Similarly, it is possible to 
find N solitary wave solutions by taking 

g (N) = j ^ e a i x + b i t + c ir (3.20) 

i=l 

where 

b t = -a * 3 (3.21) 


but the computations are very tedious for i > 3. 


4. Numerical results 

4.1. Iterative Methods Using Finite Difference Schemes 

Previously many researchers have used the finite difference methods to solve the KdV 
equation, Feng [12]. In the last decade, the CMKdV-II type equations were solved numerically 
by using split-step Fourier method [13-15]. Also parallel implementation of the split-step 
Fourier method using Fast Fourier Transform (FFT) has been studied by Taha [16] (see 
references therein). Here, in this work, the one-soliton solution of the CMKdV-II equation 
is considered. A finite interval for our numerical purposes is subjected, namely, [a,b]. The 
constants a and b can be chosen sufficiently large so that the boundaries do not affect the 
propagation of solitons. For the CMKdV-II (1.3), a numerical (finite difference) method of 
solution using iterative method is introduced, lit is approximated by using forward time 
difference scheme, U x and U xxx by the central-space difference scheme using four-points. 
Equation (1.3) becomes 


/ z «+l _ z > 

V 


11 2 

~n+ 1 , o n+1 _n+1 , n+ 1 

^m+2 T oz m+1 ^m+l m- 2 

+ 

' ~n+ 1 
^m+ 2 

9 _n+l , 9 -h+I ~n+ 1 
^m+ 1 ^m- 1 ^m-2 


12Ax 


2(Ax) 3 


where z” = z{t n ,x m ) = z{nk,mh), k = At, h = Ax. 
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0.2 



Figure 1: The modulus of the one-soliton numerical solution for the CMKdV-II equation with N = 256, 
-128 < x < 128, M = 100, 0 < t < 500. 


Multiplying both sides by 2( Ax) 3 and rearranging the terms we get 


n+l 
Zj m +2 


1 + 6(Ax) 


! (l*1 


+ z: 


n+l 
m +1 


-2 - 48(Ax) z: 


2 ^n-l 


+Z' 


+z‘ 


n+l 


2(A x) 2 

At 


+ z: 


n+l 

m -1 


2 + 48(Ax) 2 ( z "" 1 


(4.2) 


n+l 

m-2 


-1 -6(Ax) : 


(K 1 


W 7 „ 

At m 


for 3 < m < N - 2. For m = N - 1 and m = N, the backward difference scheme is chosen in 
U x and U xxx . Three more equations come from the boundary conditions, namely, L 1(a) = 0, 
L1(b) = 0, and dll/dx = 0 at x = a. Thus, N unknowns, namely, z” +1 , i = and 

N equations are obtained. Since the value of the non-linear term is known here, a system 
of linear equations is obtained. The initial guess is taken as lZ(x,0) = \/2c/a sech [y'c(x - 
Xo)]e (!00 ^ which represents a solitary wave initially at Xo moving to the right with velocity c 
and 0o is the polarization angle. The main idea is to assume that the non-linear term Iz^y 1 1“ 
is zero first and then solve the problem for whole time domain. Afterwards, this solution is 
taken and substituted for the non-linear term and solved again iteratively. The following two 
norms, namely, L^ and L 2 are used to measure the accuracy of the approximate solutions for 
stopping criteria. These norms are defined as following: 


Loo = max(||z„| - \z n \\), 

n 


l 2 = - \ z n\f, 


(4.3) 


where z„ and z n are the two consecutive new and old approximate solutions, respectively, at 
point (nAx, T) for all n, where T is the final or terminating time. FORTRAN and MATLAB 
are used to obtain the results and figures, respectively. The graph of one-soliton numerical 
solution is shown in Figure 1. 
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5. Conclusions 

In this study, Hirota's bilinear form for the complex modified Korteweg-de Vries-II 
equation is derived. One- and two-soliton solutions of the CMKdV-II equation are obtained 
analytically. One-soliton solution of the CMKdV-II equation is obtained by using finite 
difference method by implementing an iterative method. The computational cost is due to 
only finding the inverse of the matrix. The difference of two consecutive solution values 
according to the formula which is given in (4.3) is shown in Table 1 result. The convergence 
rate in the method presented above is quadratic as it can be seen in Table 1 result. It would 
be interesting to see what happens if this numerical scheme for the interaction of two-soliton 
waves for the CMKdV-II equation is applied. The numerical scheme deserves further study 
according to its application to the CMKdV-II equation. 
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